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Abstract: The number of monomers, in a monomer-dimer mean-field model with an attractive potential, 
fluctuates according to the central limit theorem when the parameters are outside the critical curve. At 
the critical point the model belongs to the same universality class of the mean-field ferromagnet. Along 
the critical curve the monomer and dimer phases coexist. 

Introduction 

Interacting particle systems described with statistical mechanics models are known to have different 
fluctuation properties on their critical points. In the mean-field ferromagnet, for instance, the sum of 
the spins centered around its mean and normalised with the square root of the total volume, converges 
toward a normal random variable (central limit theorem) away from the critical line. At the critical point 
instead a non-normal behaviour emerges, i.e. the limiting probability distribution for the sum of the spins 
centered and suitably normalised is not Gaussian [6,13]. 

In this paper we consider a mean-field system of interacting monomers and dimers where, beyond 
the hard-core interaction among particles, an attractive interaction is added to favour configurations 
where similar particles lie in neighbouring sites. The peculiar features of the presented model come from 
the combined presence of the two interactions. We show that the central limit theorem and the law of 
large numbers hold for the number of monomers at general values of the parameters. At the critical point 
instead the central limit theorem breaks down and the number of monomers centered around its mean and 
normalised with the exponent 3/4 of the total volume has a limiting density proportional to exp(— cx 4 ), 
i.e. the system exhibits the same critical behavior of the mean-field ferromagnet. We also show that along 
the critical curve the law of large numbers breaks down, due to the coexistence of the monomer and the 
dimer phases. Unlike the Curie-Weiss model, the relative weights of these phases are non-constant and 
display two contributions that correspond to the two types of interaction. 
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We provide a rigorous proof of the mentioned results by first studying the properties of the moment 
generating function for the model when the attractive interaction is zero. Here the difficulty of the problem 
stems from the fact that even in the absence of the attraction the system keeps its interacting nature and 
the equilibrium measure does not factorise. To solve this problem we use a Gaussian representation for 
the pure monomer-dimer model previously introduced in [2] which has the purpose of decoupling the hard¬ 
core interaction. When instead we consider the attractive potential we follow the Gaussian convolution 
method introduced in [7]. 

It would be interesting to further extend the results presented in this paper in the same spirit of those 
obtained for the mean-field ferromagnet in [5, 7] and also test for the same purpose other methods like 
those based on interchangeability [3,4] or those of Lee-Yang type [11]. 

The paper is organised as follows. Section 1 presents the definition of the model and the precise 
statements of the results. In Section 2 we consider the pure hard-core model and prove the law of large 
numbers and central limit theorem by using the Gaussian representation and an extended Laplace method 
(reported in the Appendix). In Section 3 we consider the hard-core model with attraction and, using the 
method of the Gaussian convolution together with the formerly introduced Gaussian representation, we 
prove the law of large numbers, the central limit theorem and their breakdown respectively along the 
critical curve and at the critical point. 


1. Definitions and Results 

Let G = (V, E ) be a finite graph with vertex set V and edge set E C {uv = {u, v} \ u, v € V, u ^ v}. 

Definition 1.1. A dimer configuration on the graph G is a set D of pairwise non-incident edges, called 
dimers: 

D C E and for each v € V there is at most one u € V such that uv £ D . (1.1) 

The associated set of dimer-free vertices, called monomers, is denoted by 

Jt G (D)~{v£V\Vu£V uv£D}. (1.2) 

We denote by the configuration space, i.e. the set of all possible dimer configurations on the graph G . 
We notice that by definition 


|^ g (D)|+2|D| = |V| VDg^ g . (1.3) 

In this paper we restrict our attention to the complete graph with vertex set V = {1,..., N} and edge 
set E = {uv | u, v € V, u v} . The corresponding configuration space will be denoted by and the set 
of monomers associated to the dimer configuration D by .#/v(D). 

A fundamental quantity is the number of monomers for a given dimer configuration D G £^j\r : 
Sn(D) := \^n(D)\. Sn(D) can be seen as a sum of N variables introducing, for a given D G 3>n 
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and for all v G V, a monomer occupancy variable 


a v (D) 


J 1, if v e ^n(D) 
0, otherwise 


Thus, one can write 


N 

S n {D) = 53 «„(£>). 

V=1 


We also define the empirical monomer density as 


m N (D) := ^S n (D ) 


(1.4) 


(1.5) 


( 1 . 6 ) 


which represents an analogous of the empirical magnetization in magnetic models. 

In [1] the authors consider a monomer-dimer model with imitative interaction on the complete graph, 
that we call Imitative Monomer-Dimer model (IMD model), defined as follows. For each integer N, the 
Hamiltonian function is 


-H n {D) := N^h-J)m N (D) + J (m N {D)) 2 ^j V£> e @ N (1.7) 

where h £ K. is the external field and J > 0 is the imitative coupling. The Hamiltonian (1.7) induces a 
Gibbs probability measure on the configurations space S>n 

» N (D) := N~\° I exp (~H N (D)) V D G ® N , (1.8) 

Zjn 

where 

Z N := N-^e~ HN(D) (1.9) 

Des>N 

is the partition function. The factor is the necessary normalisation working on the complete graph. 

As usual, the quantity 

Pn ■= ^log Z N (1.10) 

is called pressure density. 


Remark 1.2. Despite the Hamiltonian (1.7) depends only on the numbers of monomers, it is possible 
to show [1] that in our case, namely on the complete graph, any general Hamiltonian depending also 
on the number of dimers and the relative couplings is equivalent, up to a constant, to (1.7). Thus, the 
parameters h,J take into account both monomer/dimer external fields and monomer-monomer/dimer- 
dimer /monomer-dimer couplings. 

Beside the formal analogy between (1.7) and the Hamiltonian function of a Curie-Weiss model, we 
want to stress their main difference: in the former the configuration space is not a product space 
because of the hard-core constraint (1.1). 

Let us briefly recall the results obtained in [1]. We refer to the original work for the details. 
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Theorem 1.3. The thermodynamic limit of the pressure density of the IMD model is given by 

lim pn = supp(m) (1.11) 

Af->oo m 

where 

p{m) := —Jm 2 +p^°\(2m—l)J + h), (1.12) 

P ( °\h) := - I log(l - g{h)) = -1^|W - log g(h) + h , (1.13) 

9(h) e 2h + 4 — e h ) . (1.14) 

Remark l.f. We notice that, in analogy with magnetic models, one can define a general mean-field Hamil¬ 
tonian as 

-H% F (D) = Nf(m. N (D)) (1.15) 

for any bounded continuous function /. As in the case of spin mean-field models, using standard Large 

Deviations techniques, one can prove that 

lim -j- log Zff = sup (/(m) - I(m)) (1.16) 

Af->oo iV m 

where the rate function I is given by 


I(z) 


zlogz+ ±^log(l-z) + ^ -p(°)(0) lfze [0,1] 
oo otherwise. 


(1.17) 


The aim of the present work is to describe the limiting distribution of the random variable S N with 
respect to the measure /xjv, in a suitable scaling when N —> oo . From now on 5 X is the Dirac measure 
centered at x, J\f (?n, cr 2 ) denotes the normal distribution with mean m and variance a 2 and denotes 
the convergence in distribution with respect to the Gibbs measure pjy, as N — > oo. 

Let start by considering the case J = 0. The Hamiltonian (1.7) at J = 0 is a special case of the original 
problem considered by Heilmann and Lieb [9]. We introduce the following notation, 

= Z N , p^=p N , pffl = p N . (1.18) 

J =o J =0 J—0 

Setting J = 0 in Theorem 1.3 we get 


lim p$ = p^°\h) V/iG 


Thus, the pressure is analytic as a function of h and the unique value of the limiting monomer density 
is given by 

lim E ,o)(m N ) = lim p$ = p (0) (1.20) 

N—too H'N N—>o o Oil Oil 


and, using the properties of g given in [1], we get 


h >■ 


(1.21) 
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Theorem 1.5. For the IMD model at J = 0 the followings hold: 

m N 4 8 g(h) ( 1 . 22 ) 

and 

*^U*(o,>>). (1 ,3, 

We notice that, even if we are in the case J = 0, (1.23) is not a consequence of the classical central 
limit theorem, indeed Sn is not a sum of i.i.d. random variables because of the presence of the hard-core 
interaction. The proof of the theorem is in the next section. 

Let us consider now the case J > 0. It is possible to show [1] that the points where the function p 
defined in ( 1 . 12 ) reaches its maximum satisfy the following consistency equation: 

m = g((2m — 1) J + h). (1-24) 

The analysis of (1.24) allows to identify the region where there exists a unique global maximum point 
m*(h , J) of p. The resulting picture (see figure 1) is the following: the function m* is single-valued and 
continuous on the plane (h, J) with the exception of an open curve 7 defined by an implicit equation 
h = y(J). Moreover m* is smooth outside 7 union its endpoint ( h c , J c ). Instead on 7 , there are two global 
maximum points m±(J) = TOi( 7 (J), J) and m 2 (J) = 7712 ( 7 ( 0 !), J). In particular, m\ < m 2 thus they 
represents respectively the dimer and the monomer phase. The curve 7 plays a crucial physical role since 
it represents the coexistence of two different thermodynamic phases and the point ( h ci J c ) is the critical 
point of the system, whose exact value is computed in [ 1 ], 

h 

p has 1 global 
maximum point 
m* 



pf has 2 global 
maximum points 
m,<m 2 


J 

Fig. 1. The coexistence curve 7 and the critical point ( h c , Jc) in the plane (h, J). 

Outside of 7 , by differentiating the expression (1.11) with respect to the external field h, one obtains 
that the value to* maximizing p is the limit of the average monomer density to at = Sn/N with respect 
to the Gibbs measure: 

d d _ 

lim E^mAr) = lint ^tPn = 77 p(m*) = m*(h,J). 

N—too N—too Oil dh, 


(1.25) 
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We want to stress the fact that in the standard mean-field ferromagnetic model (Curie-Weiss model) the 
existence of the limiting magnetization on the coexistence curve (zero external field) is achieved by a spin 
flip symmetry argument, a property that we do not have in the present case. 

In the next sections we will prove the law of large numbers, the central limit theorem and their 
breakdowns, respectively theorems 1.6 and 1.8 below, for the distribution of Sn (suitable normalised) 
with respect to the Gibbs measure pn- 

Theorem 1.6. Consider the IMD model defined by the Hamiltonian (1.7). 

i) In the uniqueness region ( h , J) € R x R + \ 7 , we have that 

m. N . (1-26) 

ii) On the coexistence curve 7 , we have that 

'D 

niN —> Pi 5 mi + P 2 S m2 , (1-27) 

where pi = , b t = (-A;(2 - m ;))~ 1/2 and A * = ^2 p(mi), for l = 1,2. 

Remark 1.7. We notice that, on the contrary of what happens for the Curie-Weiss model, the statistical 
weights pi and p 2 on the coexistence curve are in general different, furthermore they are not simply given 
in terms of the second derivative of the variational pressure p. 



Fig. 2. We extend the definition of p\ and p 2 on a region that contains 7 , where m\ and m 2 are local maximum points 
of p, and then we compute the sign of pi — p 2 ■ The coexistence curve appears to be completely contained in the region 
Pi < P 2 - 


The first fact can be seen numerically (figure 1.7) and analytically one can compute 
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Indeed, by exploiting the formula ( p ^)" = g' = 2g(l — g )/(2 — g) (see Appendix of [1]), one can rewrite 
the ratio p\/ P 2 as 


pi _ / (2 — ?n 2 ) — 4 J m 2 (1 — to 2 ) 

P 2 V (2 — mi) — 4 J mi (1 — toi) 


(1.29) 


Furthermore, the relative weights pi have two contributions reflecting the presence of two different 


kind of interaction: the first contribution A/ is given by the second derivative of the variational pressure 
( 1 . 12 ), while the second contribution 2 — mi comes from the second derivative of the pressure of the pure 
hard-core model. 


Theorem 1.8. Consider the IMD model defined by the Hamiltonian (1.7). 
i) For ( h , J) £ (K. x M + ) \ (7 U ( h c , J c )), we have 

Sn — Nm* t> 


JV 1 / 2 


A7(V 2 ) 


where a 2 = —A 1 — (2 J) 1 > 0 and A = -£^ 2 p{m*) < 0. 
ii) At the critical point (h c , J c ), we have 


Sn — Nm c v 
N 3 / 4 > 


C ex P ( ^ x 4 ) dx 


where X c = -^ 4 p(m c ) < 0 7 m c = J c ) and C 1 = exp( 2 f x 4 )dx. 


(1.30) 


(1.31) 


2. The pure hard-core model 

A basic ingredient of all the proofs is the knowledge of the properties of the moment generating function 
of Sn w.r.t. the Gibbs measure at J = 0. However, compared with spin models, monomer-dimer models 
have an additional feature: the problem at J = 0 is itself non trivial in the sense that the Gibbs measure 
is not a product measure. We start by deriving the properties of the partition function of the model at 
J = 0 that will be used during all the proofs. 

For given u,t£l and g > 0, let us consider 

4>+^,)= E iV- |D| exp((u +] L)^(I?)) (2.1) 

-De® N 

In order to obtain an asymptotic expansion of (2.1), which allows us to obtain its scaling properties, 
we will use a connection between the monomer-dimer problem and Gaussian moments [2,14]. Following 
the same argument of [ 2 ] one finds 

Proposition 2.1. The following representation of the partition function holds 

Z n\ u +^) = /J l (*N(x)) N dx , (2.2) 

<F n (x) := (z + exp (u+J^j) exp(-y) . (2.3) 


where 
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The above Gaussian representation allows us to use Theorem A.l (see the Appendix), an extension of 
the Laplace method, to obtain a useful asymptotic expansion of Z^\u+ ^ ri ). Precisely 

Proposition 2.2. For a given u, t € K and 77 > 0 



where and g are defined respectively in (1-13) and (1.14). 

Proof. Use proposition 2.1 and check that the function Fm defined in (2.3) satisfies the hypothesis of 
Theorem A.l, with xn = e^^ u+t ^ NV ' ) g^u+t/N 11 ). By means of the stationarity condition xj v +e u+t ^ NT ’X n~ 
1 = 0, one finds log Fn(xn) = p^°\u + t/N v ) and ^ log Fn(%n) = —2 + g(u + t/N 11 ). □ 

We will show that the previous proposition gives immediately Theorem 1.5. On other hand, in the 
case J > 0 we need additional information about the convergence of to p. 

Proposition 2.3. For each k € {0,1,2,...}, Jpj4(/i) converges uniformly to ^kP^°\h) on compact 
subsets o/R. 

Proof. The location of the complex zeros h £ C of the partition function Z^\h) was described in the 
work of Heilmann and Lieb in [9]: Theorem 4.2 in [9] shows that these zeros satisfy lR(e h ) = 0, that is 
3(/i) € | +7rZ . Set U := M+i (— f, f) C C . The analytic function Z ( h ) does not vanish on the simply 
connected open set U , hence pff (h) = log Z^\h) is a well-defined analytic function on U. Moreover 
the sequence (h)) Ne N is bounded uniformly in N and in ft € K , for every K compact subset of U ; 

indeed 

|44)| < | lo S144)1 1 + ^3 

from the definition of Zfp it follows immediately 

log| 4 0 ) W| < -^log 4 } (supK(/i)) , 

and on the other hand, since zffi is a polynomial in the variable e h , using the Fundamental Theorem of 
Algebra and thank to the choice of U, it follows 



Thus, the claim is a consequence of the Vitali-Porter and Weiestrass Theorems [12]. □ 

Let now prove Theorem 1.5. For each «eR and 77 > 0 we define 

Sn, v ,u := (2.5) 

In order to prove the two statements of the Theorem 1.5, namely the law of large numbers (1.22) and 
the central limit theorem (1.23), it is enough to compute the limit of the moment generating function of 
Sn,ti,u for 77 = 1 , u = 0 and for 77 = \,u = g{h) respectively. 
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Consider the moment generating function of SV.rj.-u with respect to the Gibbs measure fi$ with 
external field h, namely for all t G M 


Des> N 


( 2 . 6 ) 


By (2.1) one can rewrite (2.6) as 


, _ -tu/N’i 44 + xn) 

VSn.v.uW ~ e - jnl,., - • 


4 ° 4 ) 

Using proposition 2.2 for the numerator and the denominator of (2.7) one gets 

44 +. 4 ) 


(2.7) 


44 ) 


exp (N(p^(h+^)-p^(h)j 


( 2 . 8 ) 


Setting 77 = 1 and u = 0 and using the Taylor expansion p(°\h+ jj) —p^(h) = jq~§^P^\h) + 0{N 2 ) 
and ( 1 . 21 ), we obtain 


Jim (j>s N ,in(i) = e t9(h) Vt G 


N—too 


(2.9) 


which implies ( 1 . 22 ). 

In the case of the central limit theorem, setting r) = | and u = g{h ), the leading order is provided by 
the Taylor expansion of p(°\h + J M ) up to the second order 

++ + 7v> = ++> + l p( ° Hh} + & + °< w + ■ 


and then we obtain 


lim 4>s , (f) = e* 2 V t G 

iV-KxC«+, g (M V 


( 2 . 10 ) 


which implies (1.23) and completes the proof. □ 


3. The model with attractive potential 

The strategy in the case J > 0 follows the general method of Ellis and Newman [5], namely, in order to 
overcome the obstacle of the quadratic term in the interaction, we consider the convolution of the Gibbs 
measure pn with a suitable Gaussian random variable. Let us start by two simple lemmas. 

Lemma 3.1. For all integer N, let Wn and Y/v be two independent random variables. Assume that 
T> 

Wn —> W, where 

Ee itw ^0 Vt G M . 

Then Yn Y if and only if Wn + Y/v W + Y . 
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Lemma 3.2. Let W ~ TV(0, (2 J) - 1 ) be a random variable independent of Sn for all TV G N. Then given 
rj > 0 and u G K, the distribution of 


W S n -Nu 
JV 1 / 2 ~v ^ jyi-v 


C N exp [nF n + u)) dx , 


where C N X = / R exp (TV Fn{-^ +u))dx and 


F n (x) := - Jx 2 + p^(2Jx + h — J) . 


Proof. Given 6 G K. , 


f W Sn ~ Nm „, 
\ TV 1 / 2-1 ? TV 1 - 1 ? 


6»| = p|-\/]VW + 5 *at G #} 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


where E = (—oo, 0TV 1_?7 + TVm]. 

The law of v^VW+Sat is given by the convolution of the Gaussian TV(0, N(2J)~ l ) with the distribution 
of S'at w.r.t. the Gibbs measure pn ■ 


NW + S N GE} = 


} 


(3.5) 


y 7 tAT 


j 


TV 


dt exp ( - i~ T t 2 ) Z^l ^ + h- J) , 


2Jt 


TV 


where the last equality follows from (2.1). Making the change of variable x = (t — Nuj/N 1 n in (3.5), we 
obtain: 

p{VNW + S n ge} =C N j dx exp ( - JTv(^ + Z ^ ^2j(^ + u) + h - J^j (3.6) 

and the integrated function can be rewritten as (3.2). □ 


The core of the problem is the convergence of the sequence of measures determined by (3.2) for suitable 
values of ry and u. Thus, we are interested in the limit of quantities of the form 

^ exp (TV F n (-^~ + ip(x) dx (3.7) 

where xf is an arbitrary bounded continuous function. Clearly, the results depend crucially on the scaling 
properties of Fn near its maximum point(s). By (3.3), (1.12) and (1.19) we know that 

lim Fn(x) = p(x), VieK. (3.8) 

N—hx 

However, the study of the asymptotic behaviour of the integral (3.7) requires stronger convergence results 
provided by propositions 2.2 and 2.3. 

Given a sequence of functions /at : ffi. —> R, for any x, y G R. we define 


Af N (x;y) := fN(x + y) - /jv(y)- 


(3.9) 
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Let p = p(h, J ) be a maximum point of p and denote by 2k the order of the first non zero derivative 
at p. Hence, making a Taylor expansion, one finds as N —» oo 

N Ap(x N~ A ; p) = -^- x ^ + o(n~^ (3.10) 

where A = &M < o- 

The next proposition relates the asymptotic behaviors of N AFm and N Ap. 

Proposition 3.3. For any x,y £ R and p > 0, 

Jim exp (^N (^Fn{xN~ v + y) — p(x N~ v + yi'j'j = c{y) (3-11) 

_1/2 

where c(y ) := (2 — g(2Jy + h — J)) . Hence, 

N(AF N (xN- r '-y)-Ap{xN- r '-y)\ -> 0. (3.12) 

V / N—><x> 

Proof. Keeping in mind the definitions (3.3), (1.12) and using Proposition 2.2 we get (3.11). Then (3.12) 
is a straightforward consequence. □ 


The next two propositions allow us to control the integral (3.7) in the large N limit. 

Proposition 3.4. Set M := max{p(a;)|a; € R}, letC be any closed (possibly unbounded) subset o/R which 
contains no global maximum points of p. Then there exists £ > 0 such that 

e~ NM J e NFN{x) dx = 0(e~ Ne ) as N ^ oo. (3.13) 

Proof. We observe that the sequence of functions (p^)jvgn is uniformly Lipschitz with constant 1, namely 
for all h, h' £ R and N € N 

IPn^W-P { N( ti )I <\h-h'i , (3.14) 

since = E^(o)(mjv) € [0,1] . From (3.14) and definition (3.3) we get 

lim svldFn(x) = —oo (3.15) 

|a?|—>-oo j\r 


and 


sup / e^^dx < oo . 

N J R 

Fixed £\ > 0, by (3.15) we can pick a number A £ R sufficiently large such that 

sup Fn(x) — M < — £i ViV G N 

X^lOa 


(3.16) 


(3.17) 


where Oa = {x € R : |or| > A}. Furthermore C \ Oa is compact (or possibly empty) and then, by 
proposition 2.3, there exist e 2 > 0 and N such that 

sup Fn{x) — M < —e 2 \/ N > N . 
xec\o A 


(3.18) 
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Thus setting e := min(£i,£ 2 ) we get 


sup Fn(x) — M < — e V TV > TV 

xGC 


(3.19) 


(3.20) 


Hence, for TV > TV, 

e-NM f e NF ^dx < e - N M e (N-l)(M-e) r e F N {x) dx 

Jc Jc 

< e~ Ne e _ ( M_£ ) [ e FN(x) dx. 

Jr 

The last is uniformly bounded in TV by (3.16) and this completes the proof. □ 

In the rest of this section d k f(x) denotes the fc th -derivative of a function / at the point x. 

Proposition 3.5. Let p be a maximum point ofp, let 2k be the order of the first non-zero derivative of 
p at p. Given 6, e > 0, there exists N e such that for all TV > TV e 

2k —1 

TV AFn{x TV“ 2 fc; p) < ^ ex^ + Lg ie x 2k \/x, \x\ < SN^ (3.21) 


3 =1 


where 


d 2k p{p) + £ su P[/i _ 5i/i+5] |<9 2fc+1 p| + £ 

^8,s : = - ,. - + 0 


(3.22) 


(2*)! (2k + 1)! 

In particular, since d 2k p(p) < 0, one can choose S, e > 0 such that Lg l£ < 0, and then the sequence of 
functions 

exp (TV AFn(x TV - ^; p)) 1 (|a;| < 5N ^) (3.23) 

turns out to be dominated by an integrable function of x. 


Proof. The Taylor expansion of Fn at the point p gives 

2k—1 


naf n(x n-* ; „ } = V a AAL N ^ x , + ** + 'LAIsM ri ,.« ,3. 

j\ (2k)\ (2k + 1)! 


24) 


3 =1 

where £ G (p, p + x N~ A). We claim that for any j G {1, ..., 2k — 1} 

d j F N (p) N^/ 2k -»• 0. 

N—too 

Indeed, by (3.12) 

n(aFn(x TVT|r ; p) — Ap(x TV~Wj p)] —>• 0 

V / N—>oo 

that is, by substituting (3.24) and (3.10) into (3.26), 


2fc—1 

E 

3 = 1 


d J F N (p) Nl _ j/2k ^ ^ d 2k F N (p) - d 2k p(p) ^ k 


hence using proposition 2.3, we get 


2/c-l 

E 

i=i 


(2fe)! 


d J F N (p) Nl _ j/2k 


O 


(tv 


N—>oo 


(3.25) 

(3.26) 

(3.27) 


N—too 


(3.28) 
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which implies (3.25) since x is arbitrary. Thus (3.25) gives the control of the terms of order up to 2k — 1 
in (3.24). The last two terms in (3.24) can be grouped together observing that Ixl 2 *" 1 " 1 < x 2 k SN A; then 
the estimate (3.21) is obtained using the uniform convergence of d 2 k F N , d 2 k+ 1 F N on the compact set 
[p — 6 , /i + 5], which is guaranteed by proposition 2.3. □ 


Let now prove Theorem 1.6. We denote by A4 = {mi}i=i,...,p the set global maximum points of p 
and let ki and A i be as in (3.10). Set M := ma x m p(rn) = p(pi) for each l = 1,..., P. From the analysis of 
p and using the properties of the function g (see [ 1 ]), it turns out that ki do not depend on l and precisely 


(M,k) = { 


({m*(/i, J)}, l) 
({m c }, 2 ) 


if (h, J) G (R x M+) \ (7 U (h c , J c )) 
if (h, J) = (h c , J c ) 


^({toi(J),to 2 (J)},1) if (h, J) G 7. 


(3.29) 


The argument described below applies in all the cases proving respectively (1.26) and (1.27). Keeping 
in mind (3.29), we proceed with the computation of the limiting distribution of the monomer density 
m.N = Sn/N. By lemmas 3.1 and 3.2 with 77 = 0 and u = 0, it suffices to prove that for any bounded 
continuous function ip 




NF n (x 


] d 


x 


p 


E 

1=1 


E b, 


1 =1 


(3.30) 


For each / = 1,..., P let 61 > 0 such that the sequence of functions (3.23), with /i; in place of p, is 
dominated by an integrable function. We choose 6 = min{4/ | l = 1,. .., P}, decreasing it (if necessary) 
to assure that 0 < 5 < min{|/j; — p s \ : 1 < l 7 ^ s < P}. Denote by C the closed set 

p 

C := R\ {J(vi - S,ni + 5) ; 

1=1 


by proposition 3.4 there exists £ > 0 such that as N —> oo 


—NM 


J e NFN ^ip{x)dx = 0 (e~ Ne ) . 


(3.31) 


For each ! = l,...,Pwe have 

r/ii+8 


N™ e ~ NM / e NFN{x) ip{x) dx = 

J ui—6 


' IH—S 
= e N(F N (vn)-M) 


i exp (nAFjst(wN ss ; pi M ip(wN 2 * + pp) dw 


(3.32) 


J |n;|<<5iV 

where the equality follows from the change of variable x = pi + wN~^k and AFn is defined in (3.9). 
Since M = p{pi ), from (3.11) we know that 


1 


1 


lim p N(Fn(im)~M) = _ 

- g(2Jpi +h-J) V 2 ^ di 
where the last equality follows from the fact that pi must satisfy the consistency equation (1.24). 


(3.33) 
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By Proposition 3.5 we can apply the dominated convergence theorem to the integral on the r.h.s. of 
(3.32), then by (3.12) and (3.10) we obtain 

r^i+s 

N * e~ NM I 

N—>oo 


rl±l+5 

lira N™ e~ NM / e NFN{x) ip(x)dx = 

hi-s 




l = f exp ( 

~ IH J R V 


^ l ,2 k 


(3.34) 


(2fc)! 


ur k ) ip(ni)dw . 


Making the change of variable x = w(—Xi) 2k in the r.h.s. of (3.34) and using (3.31) we obtain 


lim N 2k e 

N—>oo 


^ -NM 


0 N F n (x 


P i - 

] ip{x)dx = /2 _ (-A/ exp ^ 

^^ * r'l 




(2fc)! 


dx . 


(3.35) 


The analogous limit for the denominator of (3.30) follows from (3.35) by choosing ip = 1. This concludes 
the proof of the Theorem 1.6. □ 

Let now prove the Theorem 1.8. Keeping in mind (3.29), let us start by proving the following 


/ R <^(pj, *”)«*><** 


/ exp ( N Fn (x N 2 k +m*))dx [ exp ( \ x 2k ^\ dx 

Jr ' Jr '( 2 k )' ' 


(3.36) 


for any bounded continuous function ip. We pick 5 > 0 such that the sequence of functions (3.23) is 
dominated by a integrable function. By proposition 3.4 there exists e > 0 such that as IV —► oo 


-NM 


[ i exp (nFn(xN 2k + to*)') ip(x) dx = 0(N 2k e Ne ) (3.37) 

J\x\>8nJh ' ' 


where M = rna x m p{m). On the other hand as \x\ < SN 1 /^ 

e ~ NM f i exp (nFjv^xN-J* + m*)') ip(x) dx = 
J\x\<8nJh ' ' 


• | a? | <<5iV 2 
= e (Fi»(m*)-M) 


i exp ^ N AFn[ xN 2k -,m*)^ ip(x) dx . 


(3.38) 


lim e~ NM 

N—too 


' \x\<. 6 N 2 ic 


exp ^ N Fn{xN 2k + m*)^j ip(x) dx 


J\x\<8N2H 

Thus, by proposition 3.5 we can apply the dominated convergence theorem, and then by (3.33), (3.12) 
and (3.10) we obtain 

(3.39) 

which, combined with (3.37), implies (3.36). 

For k = 2, by lemmas 3.1 and 3.2 with 77 = 1/4 and u = m *, the convergence (3.36) is enough to 
obtain (1.31). 

For k = 1, by lemmas 3.1 and 3.2 with 77 = 1/2 and u = to*, since W ~ A/"(0, (2J) _1 ), the equation 
(3.36) implies that the random variable Sn converges to a Gaussian whose variance is a 2 = (—A )^ 1 — 
( 2 J) _1 , provided that 

\ 4- 9 7 

(3.40) 


(-A )- 1 - ( 2 J )- 1 = > 0 


-2AJ 

where A = -§^p{m*). Considering the function g defined in (1.14), we have that J^p(m*) + 2J = 
(2J ) 2 g’(2Jm* + h — J). Since g’ > 0 and A < 0 the inequality (3.40) holds true. □ 
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A. Extended Laplace Method 


The usual Laplace method deals with integrals of the form 


(■ ip(x)) n dx 


as n goes to infinity. In this appendix we prove a slight extension of the previous method where ip can 
depend on n. Other results in this direction can be found in [8,10]. 


Theorem A.l. For all ?j€N let ip n : K — > R . Suppose there exists a compact interval A'Cl such that 
ip n > 0 on K, so that 

ip n (x) = \/x £ K . 

Suppose that f n £ C 2 (K) and 

a) f n — > f uniformly on K ; 

n—>o o 

b) /" —> f" uniformly on I\ . 

n—>oo 

Moreover suppose that: 


1) maxjf /„ is attained in a point x n £ int(K); 

2) limsup,^^ (sup R \ A - log \ip n \ - max A /„) < 0 ; 

3) maxjf / is attained in a unique point x £ K; 

4) f"{x) < 0 ; 

5) limsup n _ > , 00 f R \ip n (x)\ dx < oo . 


Then, 


..(a)) "da ~ e nfn(x„) 


27T 


In the proof we use the following elementary fact: 


(A.l) 


Lemma A. 2. Let (f n ) n be a sequence of continuous functions uniformly convergent to f on a compact 
set K . Let ( I n ) n and I be subsets of I\ such that madist(a:, y) —> 0 as n —> oo . Then 

• max/ n f n ->• max/ / ; 

n—> oo 

• argmax 7 f n —> argmaxj / , provided that f has a unique global maximum point on I. 

71 n—)• oo 
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Proof (of the Theorem A.l). Since x n is an internal maximum point for f n (hypothesis 1), f' n {x n ) = 0 
and for all x G K 

fn{x) = fn{x n ) + ^ fn(fx,n) (X - X n f with G (x n , x) C K . (A.2) 

Fix £ > 0. Since /" —► f" uniformly on K, there exists N e such that 

n—>■ oo 

\fn(0 ~f "(01 < £ V^GA' Vn > N s . (A.3) 

Since f" is continuous in x, there exists S e > 0 such that B(x,S e ) C I\ and 

\f"{0-f"{x)\<e VH:\d-x\<8 e . (A.4) 


By the lemma A.2 x n —> x, because x is the unique maximum point of / on K (hypothesis 3). Thus 

n—to o 

there exists N$ e such that 

\x n -x\ < y \/N>N Se . (A.5) 

Therefore for n> N e V N$ e and x G B(x, 5 e ) it holds: 

(A.5) S £ 

|£:r,n %\ — x| ~b |x x\ ^ \x n x\ H - \x x\ <C ^ = ^ > 

l/A'fo.n) - /"(®)| < |/"(£*,n) - /"(&,„)! + |/"(£*,n) - f"(x) I (A ' 3) < (A - 4) £ + £ = 2 £ . 

By substituting into (A.2) we obtain that for n> N e V N$ e and x G B(x,S e ) 

< fn(x n ) + \ (f"(x) + 2e) (x - X n ) 2 

(A.6) 

> fn{x n ) + \ (f"(x) ~ 2e) {x - Xn ) 2 
Now split the integral into two parts: 

[ (ipn(x)) n dx = [ e n f n(x„) dx + f (^ n (x)) n dx . (A.7) 

Jr J B{xnAe) JM\B(x n ,5 £ ) 

• To control the second integral on the r.h.s. of (A.7) we claim that there exists rj$ e > 0 and Nf such 

that 

log \ip n (x)\ < f n {x n ) — Vs e V x G R \ B{x n , 8 e ) \/N>Ng e ; (A.8) 

namely limsup^^ sup xeR \ B(4n;5£) log \if> n (x) \ - f n (x n ) < 0 . Indeed: 


fn(x) 


limsup sup log \4>n(x)\ - fn(x n ) = 


n-¥ oo i6l\B(i„,j £ ) 


limsup sup f n (x) - fn{x n ) V 

y n-¥ oo x£K\B(xn,S e ) ) 


( limsup sup log|^ n (x)| 

n-¥ oo xeM\K 



sup f(x) - f(x) 
ixGK\B(x,S e ) , 


V 


( limsup sup log|^ n (x)| 

n->oo x£R\K 


fn(Xn) 


where the last identity holds true by the lemma A.2. 

Moreover sup^g^-^^ ^) f(x) — f(x) < 0 since x is the unique maximum point of the continuous function 
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/ on the compact set K (hypothesis 3); while limsupj^^ sup xeR y /f log \ip n (x)\ — f n {xn) < 0 by the 
hypothesis 2. This proves the claim. 

Now using (A.8) and the hypothesis 5, there exist C and N such that for all n>iVV N$ 


IR\B(x n ,S e ) 


|^> n (ar)|"dx < e^ n TIAT™) vs e ) f \ip n ( x )\dx 

Jr 


(A.9) 


• To study the first integral on the r.h.s. of (A.7), choose e € (0, £o], where f"{x) + 2e 0 < 0 (hypothesis 
4)- By (A.6), since we can compute Gaussian integrals, we find an upper bound: 


^/nfzOdx < e «/n(£n) / gf (/"(*)+2e) (*-i„) 2 fa 


'B(x n ,6 e ) 


_ e "/n(*n) 


_ e nf„(x „) 


V-§(/"(*)+ 2e) 


e - * dx 


27T 


—n (f"(x) + 2e) 


(A.10) 


and a lower bound: 


/ e n/„(x) d x > e nfn(x„) / e f(/"(x)+2e)(x-i n ) dx 

I B(x n ,S £ ) J B(x n ,5 £ ) 


= e nMi 


— e nfn(Xn 


1 f 

\j~ f {f"{x) + 2e) 7s(o,<5 £ ^/-f (/"(®)-2e)) 
/ 27T 


-n(f"(x) - 2e) 


(1 + UJ n ,s,S e ) 


(A.ll) 


where Wn,e,< 5 e —> 0 as n —» oo and £ is fixed. 

In conclusion, by (A.7), (A.9), (A.10), (A.ll) we obtain that for £ G (0,£o] and n > N e VN$ e VNVN^ 
it holds: 


Jr (V’n(z)r dx 




-nf"(x) 


< 


/"(£) 


f"{x) + 2e 


2tt 


/ /"(*) 

n->oo V /"(x) + 2£ e->0 


1 


and: 


/r (V’n( a; )) n dx 


lfn{x 


> 


/"(*) n , ^ r .l n f"(x) r -r 

(1 “I” ^n,£,6 £ ) L/ y ~ ^ 


-nf"{x) 


f"{x) - 2£ 


27T 


/ ///( ^ 

ra->oo V f"(x) — 2e e->0 


1 


hence (A.l) is proved. □ 
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